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<N ; Abstract 

We study the four-body r — > v T l I X^ decays where I = e or u and X = tt, K, p and K* 



mesons. These decay processes violate the total lepton number (|AL| = 2 ) and can be induced by 
the exchange of Majorana neutrinos. We consider an scenario where these decays are dominated by 



the exchange of only one heavy neutrino which produces an enhancement of the decay amplitude 

-L ' via the resonant mechanism. Searches for these novel decay channels with branching fractions 

sensitivities of O(10~ 7 ) can provide constraints on the parameter space of the Majorana neutrinos 



which are stronger than the ones obtained from AL = 2 decays of charged pseudoscalar mesons. 
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I. INTRODUCTION 



Current evidence from oscillations experiments [lj allows to conclude that the involved 
neutrinos are very light massive particles and that their flavors are mixed. One of the most 
intriguing and still unsolved questions in particle physics is to elucidate if neutrinos are 
Dirac or Majorana fermions [2j. The Majorana nature of neutrinos can be established in 
the simplest way through the observation of processes where the total lepton-number L is 
violated by two units QAL| = 2), a property that emerges from the non-invariance of the 
neutrino mass term |3j-l7| under global phase transformations of Majorana fields. 

Up to now, most experimental efforts have focused on searches of neutrinoless double beta 
nuclear decays [8L which is by far the most sensitive AL = 2 channel. The non-observation 
of these decays 9Ml3l| has provided very strong constraints on the existence of very light 
Majorana neutrinos, and has established direct upper bounds on the effective Majorana mass 
of the electron-neutrino (m ee ) at the sub-eV level 13j (here we define (mw) = J2i UuUiumi, 
where 1,1' = e,fi,r and i = 1,2,3 labels the neutrino mass eigenstates; see Section II for 
notation). Direct bounds on other entries of the effective Majorana mass matrix are very 
ooorly known [13J, but indirect upper limits can be obtained by combining oscillation data 
14| . cosmological bounds |l5Ml8| and tritium beta decay [191]. In turn, these indirect bounds 
on (mui) can be used to predict other AL = 2 decays, such as same-sign dileptons produced 
in r lepton or K, D and B meson decays. The predicted rates in this light Majorana neutrino 
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221 ] and beyond the sensitivities of current 



scenario turn out to be extremely suppressed 
and future superflavor factories. 

A nice explanation for the very light scale of neutrinos can be found in the existence of 
additional heavy right-handed neutrinos via the see-saw mechanism 23j-|29|. As a remnant 
of lepton-number violating Majorana mass terms, the couplings of lepton charged currents 
can induce AL = 2 processes when expressed in the basis of Majorana neutrino mass 
eigenstates. However, the exchange of very light or very heavy Majorana neutrinos in 
these decays are strongly suppressed and usually also lead to unobservable rates 
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22]. 



An alternative scenario is provided by the so-called resonant mechanism [30] which can 
produce large enhancements of the AL = 2 transition amplitudes if the masses of exchanged 
Majorana neutrinos are accessible to the energy scales of the physical processes. In this case, 
the non-observation of lepton number violating decays can be turned out into significant 



constraints on the mixings and masses of Majorana neutrinos. Let us note that such heavy 
neutrinos in the range of a few keV to a few GeV can play an important role in cosmology 
and astrophysical processes |3l|, without conflicting neutrino oscillation data. For instance, 
some extensions of the Standard Model incorporating right-handed singlet neutrinos provide 
a good candidate for dark matter in the form of a stable sterile neutrino in the range of 
a few keVs 31]; also, such models contain additional heavier sterile neutrinos with masses 
of 0(1 GeV) which can explain the baryon asymmetry of the universe 3jJ. Alternatively, 
Majorana neutrinos with masses in the range of a few hundred MeV to a few hundred 



GeV can be generated dynamica 
electroweak symmetry breaking 



ly in extended technicolor model realizations of dynamical 



32|. 



At low energies, the resonant enhancement scenario has been studied in several decays 
of pseudoscalar mesons and tau leptons. The three-body decays t~ — y l + Mj~ M^ . K + — y 



n-l + l' + , D + -> M~l + l' + , B + -> M~l + l ,+ have been considered in Refs. |20j, |30|, Mw\ 
(charged conjugated modes are implied in all channels). In the case of the r lepton decays 
some constraints can be derived on the product of two different mixing angles (for instance 
IV^jvKvI) as a function of the neutrino mass rriN |30|, while the decays of pseudoscalar 
charged mesons allow to constrain also the individual mixing angles, for instance \Vin\, as a 
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function of the neutrino mass 

of the four body decays B° — y D~1 + 1 + tt 



36[. Very recently, we have reported the first calculation 
which are expected to provide complementary 



constraints to the three-body decays of their charged counterparts. At higher en ergies, 
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the 



46[ and in top quark decays [38 
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48| 



production of same-sign dileptons at colliders [30 
has also been considered in the literature. 

Searches for AL = 2 three-body decays have been carried out by several experiments and 
the upper limits on the branching fractions can be found in Refs. 49j-[52j]. New upper limits 
on branching fractions of lepton number viola ting decays of charged B mesons have been 
reported recently: (1) the Belle Collaboration 53J has obtained B(B~ —y l~l~D + ) < 10 -6 
(I = e, ft) at the 95% confidence level (CL); (2) the BABAR Collaboration 54J has reported 
upper limits at the 90% CL for B + —y h~l + l + (h = tt/K, I = e//i) of the order of a few 
times 10~ 8 ; (3) the LHCb Collaboration has obtained results on the B~ —y X + fi~n~ decays 
(X = D, D*, D s , 7T, K) with upper bounds ranging from 10 -6 to 10~ 8 [5JJ |56]; in addition, 
an upper limit has been reported for the four-body decay B(B~ — y D°tt + fi~ fi~) < 1.5 x 10 -6 
at the 95% CL [56(. Finally, searches for AL = 2 decays of the r lepton have been reported 



by the Belle collaboration on six different r~ — > l + M^ M^ decay channels, with upper limits 
on branching ratios of the order of 10 -8 [57j. 

In the present paper we study the AL = 2 tau lepton decays r 1 * 1 — > z/ T / ± / ± X =F (with 
I = e or fi, and X = it, K, p or K* meson) within the scenario provided by the resonant 
Majorana mechanism. These decays allow to derive bounds on the |Vzjv| (/ = e /i) mixings, 
contrary to the case of three-body r lepton decays which only allow to derive bounds on the 
product |V/jvK-iv|- Given the clean environment provided by r lepton decays, these bounds 
on the Majorana neutrino mixings are free from hadronic uncertainties that are intrinsic 
to decays of pseudoscalar mesons. Therefore, these novel decay channels allow to derive 
constraints on the mixings that are complementary to those obtained from tau lepton and 
meson decays. The large sample of r lepton pairs (~ 10 10 ) that are expected to be recorded 



at the superflavor factories 58|, |59j, makes very attractive the study of these lepton number 
violating processes. 



II. CHARGED CURRENTS OF MAJORANA NEUTRINOS AND KINEMATICS 

The Feynman diagram corresponding to the T~(p) — > v T (pi)l~ (p2)l~ (j>z)X + (p 4 ) decays, 
where X can be a pseudoscalar or a vector meson, is shown in Figure 1. The letters within 
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we 



brackets label the momenta of each particle. Following the definitions given in Ref. 

can write the differential decay rate (in the rest frame of the decaying particle of mass M ) 

in terms of five independent kinematical variables (see conventions in Figure 2): 



*/3l2# 
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dT = — — — — - |7W I 2 • — -dsudsud cos 6\d cos Ozdd) , (1) 

4(47r) 6 M 3 ' ' n\ 

where Si 2 = {p\ + p 2 ) 2 and s 34 = (p 3 + p 4 ) 2 denote the invariant masses of the 12 and 34 

particles, while (#i, 6*3, <p) are angular variables defined in Figure 2 38|,|60|. The n! factor 

in the denominator of Eq. (1) accounts for identical particles in the phase space, \Ai\ 2 is the 

spin-averaged and properly antisymmetrized (under exchange of identical leptons) squared 

amplitude, /3 12 (^34) is the velocity of particle 1 (particle 3) in the center of mass frame of 

particles 1 and 2 (3 and 4) and X = \(p 2 — sp . — s 34 ) 2 — 4S12S34] 1 / 2 . 



Similarly to previous studies 
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471 ]. we add n right-handed singlets N bR (b = 
1,2, •••n) fields to the usual three left-handed SU(2) lepton doublets L T aL = (u a , Io)l, 
(a = 1,2,3) of the Standard Model. In terms of the neutrino mass eigenstates obtained 




FIG. 1. Feynman diagram for the r 
Major ana neutrino is denoted by v^ 



v T l I X + decay, where X = n, p, K or K*. The 



from the diagonalization of the Dirac and Majorana mass terms, we can write the charged 



current interactions of leptons as follows 

r 3 



30: 



c? = -^wi 



J2 J2 U lmPm-fPLl + J2Y1 VlmN^PJ I + h. C. 



(2) 



l=e m=l l=e m=l 

where Pl = (1 — 7s)/2 is the left-handed chirality operator, g is the SU(2)i gauge coupling, 
tj) c = C%j) T is the charge conjugated spinor, and U[ m (VJ m ) denotes the mixings of light 
(heavy) neutrinos; the subscript m refers to the basis of mass eigenstates obtained from the 
diagonalized Majorana mass term for neutrinos 30]: 






~2 [ S m ™ V mhV C mR + Y^ m m N mL N mR J + h. C. . 

\m=l m=l / 



(3) 



As in previous studies, we will assume that only one heavy neutrino with mass win and 
charged current couplings Vin to leptons, dominates the decay amplitudes via the resonant 
enhancement mechanism. This scenario is useful to simplify the analysis of the parameter 
space, and it can be accomplished if the spectra of heavy neutrinos is such that only one of 
them falls in the mass region that is accessible in the decay under consideration (mi + m^ < 
mjv < m T — mi in the present case). 

III. LEPTON NUMBER VIOLATION IN FOUR-BODY r LEPTON DECAYS 

Following the convention of momenta defined in the previous section, we can write the 
(properly antisymmetrized) decay amplitude for the AL = 2 decays of the r lepton as 




FIG. 2. Kinematics of a generic four-body decay in the rest frame of the decaying particle, 
J2i=i Pi = ®- We have defined pij = Pi + Pj, such that P12 + P34 = 0. 

follows: 

M = (^ F V^ N m N u Vr {p 1 )'fP L u T {p) ■ u(p 2 ) [V N (p 2 h lllu + V N (p 3 h ulll ] u c (p 3 ) (V u ) . (4) 

Here, Gp denotes de Fermi coupling constant, V uq (with q = d or s) is the entry of the 
Cabibbo-Kobayashi-Maskawa matrix for the hadronic vertex, Vin is the neutrino mixing 
defined in Eq. (2) and m^ denotes the mass of the heavy neutrino. In the case of two 
different charged lepton flavors in the final state, which we do not consider in this paper, 
we should replace Vf N — > VJ=atV^at. As it was stated before, we consider that only one 
heavy neutrino iV dominates the decay amplitude. The Lorentz-vector in Eq. (4) becomes 
V u = ijpp\ when X is a pseudoscalar meson and V u = fv m v eU {p±) when X is a vector 
meson. In our numerical evaluations, we will use the following values of the meson decay 
constants (all given in MeV units): /„- = 130.4, f% = 156.1 from Ref. [49_j, and f p = 216 
■^eV , fx* — 205.4 where obtained from the measured rates of r — > Vv T decays quoted in 
491 ] . The lifetime of the r lepton and the values of the quark mixing angles were taken also 
from [491 ]. 

In the expression for the decay amplitude we have introduced the factor 

1 



V N {pi) 



{Q - Pi) 2 - m N + im N r 



N 



(5) 



where Q = p — pi = P2 + Pz + Pi- In this expression, Tn represents the decay width of the 
heavy neutrino; it allows to keep finite the amplitude when the heavy neutrino is produced 
on-shell, (Q — Pi) 2 = m 2 N . For a given mass m^ of the heavy neutrino, its decay width can 
be obtained by adding up the contributions of all its decay channels that can be opened at 
the mass m N |30| : 

IV = ^T(iV -+ f)9(m N -£>,<) > ( 6 ) 

/ i 

where m^ in the argument of the step function are the masses of the final state particles 
in the neutrino decay channel /. The dominant decay modes of the neutrino in the range 
of masses that are relevant for the resonant r lepton decays are induced by the exchange of 
W ± , see Eq. (2), and Z° gauge bosons: iV -)■ FP ± , u t P°, PV ± , u t V , lflfu h , u h l^l+, and 
v\ x vi>, where I, l\, li = e, /x, and P (V) denotes a pseudoscalar (vector) meson state. The 
expressions for the partial decay rates of these channels can be found in Appendix C of Ref. 



30] 



As it can be checked from Figure 4 in Ref. 38), the decay width Tn varies between 10 -20 
GeV and (at most) 10~ 14 GeV for values of neutrino masses that are relevant for resonant 
t~ — > v T l~l~X + decays. These numerical values are indeed upper limits and were obtained 
by assuming the bounds on the mixings of the heavy neutrino with the three charged leptons 



as reported in Ref. 61| . namely: 

\V eN \ 2 < 3 x KT 3 , |V>| 2 < 3 x 10~ 3 , \V tN \ 2 < 6 x KT 3 . (7) 

In other words, the neutrino decay width is so tiny that, for our purposes, we can use the 
narrow width approximation, 

lim V N (pi) = -mS ((Q - Pi) 2 - m 2 N ) , (8) 

to convert the five-dimensional integral in Eq. (1) into a four-dimentional one. The branch- 
ing ratios are then obtained by using the Montecarlo code VEGAS to perform numerically 
the four- dimensional integration. 

So far, no experimental searches have been reported for the t^ 1 — > z/ r / ± Z' ± A =F decays. 
With the large data sample of r lepton pairs that are expected at superflavor factories 
59| | . we may expect that sensitivities at or below the 10~ 7 level may be easily reached 
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for the branching ratios of these decay channels. Just to illustrate the potencial of r lepton 
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FIG. 3. Branching fractions for r — > u T e e X^ decays as a function of mj^ 
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FIG. 4. Same description as in Figure 3 for di-muon channels. 

decays to constrain the parameter space of the heavy neutrino, in Figure 3 (Figure 4) we 
show the calculated branching ratios for the di-electron (respectively, di-muon) channels as 
a function of the Majorana neutrino mass m^ by using the upper bounds shown in Eq. (7). 

When upper bounds on the branching ratios of r — > v T Xll decays become available, we 
will will be able to get constraints in the \Vi^\ vs. m^ plane, as done for example in Ref. 
[301 ] . These contraints are derived by noticing that, in the narrow width approximation, the 
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FIG. 5. Exclusion regions in the |V^»jv|-mAr plane, by assuming upper bounds on B(r 
v T e~e~X + ) decays of order 10~ 7 . 

dependence of the branching ratios upon the mixing angles is as follows: 



Bfr* -> v T rrx T ) 



IT/ I 2 



'IN\ 



■h 



(9) 



h\V eN \ 2 + f 2 K N \ 2 

where fi (i = 1, 2, 3) depends upon the relevant coupling constants, phase-space integrals 
and the mass of the heavy Majorana neutrino. In order to illustrate the constraints that 
can be gotten from the experimental searches, we will assume upper limits of O(10 -7 ) for 
the branching ratios of different decay channels and we set |V^v| = {V^nI 62|- In Figures 5 
and 6 we plot the exclusion regions (region above the plotted curves) in the |V/jv| 2 vs. m N 
plane for the di-electronic and di-muonic channels, respectively. 

As we can observe from Figures 5 and 6, both leptonic decay channels offer different sen- 
sitivities to the parameter space of the additional heavy neutrino. The searches of r decays 
with different pseudoscalar and vector mesons in the final state would allow to constrain 
a larger and complementary region of the parameter space, although for the assumed 10~ 7 
upper limits most of the excluded regions can be obtained from the strangeness-conserving 
channels. On the other hand, a comparison of our results with Figures 9 and 11 of Ref. 30 1 
shows that the r lepton decays considered in this paper can provide stronger constraints on 
the \Vin\ vs. itln parameter space than the ones coming from decays of D and B meson 
decays, at least in the region of neutrino masses where these decays overlap. The results 
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FIG. 6. Exclusion regions in the |V^7v|-mAr plane, by assuming upper bounds on 5(r — > 
is T fi~fi~X + ) decays of order 10 . 

of the present paper indicate that AL = 2 decays of the r lepton can provide competitive 
constraints not only for the the product iVz/yV^/yl of neutrino mixing angles, as is the case 
of three-body decays, but also over | Vw\ without further theoretical uncertainties related to 
hadronic form factors or loop effects 21]. 



IV. CONCLUSIONS 



In this paper we have studied the four-body r ± — ¥ v T l V^X^ decays, where I = e or p 
and X = 7i~, p, K and K* mesons. We consider a framework where the heavy Majorana 
neutrino that mediates these AL = 2 decays can enhance the decay amplitudes via the 
resonant mechanism and that the contribution of only one heavy neutrino dominates the 
decay amplitude. 

We have found that these novel four-body decays, together with the three-body r ± — > 
l^M^M^ decays previously studied by other authors, can provide a more complete set of 
constraints on the parameter space associated to the mass and mixings of the Majorana 
neutrino. One important advantage of these four-body r lepton decays is that they are free 
from the hadronic uncertainties associated to the decays of pseudoscalar charged mesons and 
depend only on well known decay constants of pseudocalar and vector mesons. By assuming 

10 



experimental sensitivities of O(10~ 7 ) for branching ratios of different channels at superflavor 
factories, we find that the AL = 2 four-body decays of r leptons can provide constraints 
for mixing angles |V/tv| 2 ~ 10~ 3 to 10~ 4 , which are similar or better that the ones obtained 
from B and D meson decays. 
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